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We investigate linear and nonlinear transport across single-walled carbon nanotube quantum 
dots weakly coupled to spin-polarized leads. We consider metallic tubes of finite length and small 
diameter, where not only forward scattering contributions of the Coulomb potential, but also short- 
ranged processes play an important role. In particular, they induce exchange effects leading for 
electron fillings 4n+2 either to a non-degenerate groundstate of spin 5* = or to a triplet groundstate. 
In the linear regime we present analytical results for the conductance - for both the S — and 
the triplet groundstate - and demonstrate that an external magnetic field is crucial to reveal the 
spin nature of the groundstates. In the nonlinear regime we show stability diagrams that clearly 
distinguish between the different groundstates. We observe a negative differential conductance 
(NDC) effect in the S — ^ groundstate for antiparallel lead magnetization. In presence of an 
external magnetic field, spin blockade effects can be detected, again leading to NDC effects for both 
groundstates. 

PACS numbers: 73.63.Fg, 72.25.-b, 73.23.Hk, 85.75.-d 



INTRODUCTION 



Since their discovery by S. lijima and T. Ichihashi^ in 
1993, single-walled carbon nanotubes (SWNTs) have at- 
tracted attention due to their remarkable electronic and 
mechanical properties^^^. At low energies, they repre- 
sent an almost perfect realization of a one-dimensional 
(ID) system of interacting electrons with an additional 
orbital degree of freedom due to the sublattice structure 
of graphene. Accounting for spin and orbital degrees of 
freedom implies that for nanotubes a shell structure is 
expected, where each shell can accommodate up to four 
electrons. In the absence of Coulomb interaction the en- 
ergy levels are spin degenerate, while the orbital degen- 
eracy is usually lifted due to the nanotube finite length. 
Coulomb interactions, however, modify this picture. The 
sublattice structure of graphene gives rise to a distinction 
between electron interactions on the same and on differ- 
ent sublattices. Therefore, besides the long-ranged for- 
ward scattering processes, also short-ranged interaction 
processes play a role in small diameter tubes^^. These 
short-ranged interactions cause in finite size nanotubes 
exchange effects leading for a tube filling of 4n + 2 to a 
groundstate with either total spin = or 6* = ^ (a 
triplet)^. Signatures of the exchange interactions have 
indeed been inferred from stability diagrams of carbon- 
nanotube-based quantum dots^^^. In particular it was 
shown by Moriyama et al.^ that an applied magnetic field 
can be used to reversibly change the groundstate from the 
singlet to one of the triplet states. 

Recently, carbon nanotubes have also attracted much 
attention for their potential applications in spintronic 
devices^^. They are particularly interesting because they 
have a long spin lifetime and can be contacted with fer- 
romagnetic materials. Indeed, spin-dependent transport 
in carbon nanotube spin valves has been demonstrated 
by various experimental group^~— , ranging from the 
Fabry-Perot^^^^^ to the Kondo regime^. 
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FIG. 1: Single-electron-tunneling setup of a single- walled car- 
bon nanotube (SWNT) which is weakly coupled to source and 
drain contacts. The contact magnetization may either be par- 
allel or antiparallel as indicated by the arrows. The gate elec- 
trode allows to shift the chemical potential inside the SWNT. 



From the theoretical point of view, spin-dependent 
transport in interacting SWNTs has been discussed so 
far in the limit of very long nanotubes^^, for tubes in 
the Fabry-Perot regime^^ and for SWNT-based quantum 
dotsi2r— . In the three latter works the characteristic 
four-electron shell-filling could be observed in the sta- 
bility diagrams. In^^ however, focus was on medium- 
to-large diameter SWNTs where exchange effects can be 
neglected. The studies in^^-^^ are based on the theory 
by Oreg et al.^, where exchange interactions are treated 
on a mean-field level, and focus predominantly on shot 
noise^^ and cotunneling^^ effects. 

In this work we generalize the previous investigations of 
Ref to include the short range Coulomb interactions 
causing exchange splittings of the six otherwise degen- 
erate (at vanishing orbital mismatch) 4n + 2 - filling 
groundstates. The leads are either parallel or antipar- 
allel spin-polarized and weakly coupled to the SWNT, 
see Fig. [TJ In the low bias regime we derive analytical 
formulas for the conductance for both large and small 
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orbital mismatch corresponding to an 5* = and S = H 
groundstate, respectively, at 4n + 2 filling. In the high 
bias regime we numerically calculate the stability dia- 
grams for the two possible groundstates. We show several 
differences in transport between parallel and antiparallel 
lead magnetization, as e.g. a negative differential conduc- 
tance (NDC) effect occurring only for the S = ground- 
state and antiparallel magnetization. We further include 
in the calculations a parallel magnetic field leading to a 
Zeeman splitting for all states with total spin unequal to 
zero. It is then possible to observe spin blocking effects 
due to transport channels that trap the system in the 
triplet state with Sz = —fi- Performing a magnetic field 
sweep, a groundstate change may be obtained as it has 
been shown experimentally^. 

The paper is organized as follows. In section II we dis- 
cuss the relevant features of the low energy Hamiltonian 
of interacting SWNTs with special focus on the filling 
4n + 2. In section III we describe the set-up and method 
used to study spin-dependent transport in the sequential 
tunneling regime. Finally, in section IV, we present our 
results for the conductance, while in section V we focus 
on the nonlinear (finite bias) regime. 
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FIG. 2: The dispersion relation of a noninteracting SWNT 
with open boundary conditions. It is characterized by two 
linear branches, r = =b, of slope ±hvF determined by the 
Fermi velocity vf- The quantities eo and ca are the intraband 
level spacing and the orbital mismatch energy, respectively. 

Fermi-sea, correspondingly reads 



II. THE INTERACTING LOW ENERGY 
SPECTRUM 

A. The interacting Hamiltonian 

The starting point for a microscopic, but still analyti- 
cal, treatment of SWNTs is a tight-binding ansatz for the 
wavefunction of the 2pz - electrons on the graphene hon- 
eycomb lattice. Including nearest neighbor hopping ma- 
trix elements it yields an electron-hole symmetric band- 
structure with a fully occupied valence band and an 
empty conduction band. Since the two bands touch at 
the cornerpoints of the 1st Brillouin zone, the Fermi- 
points, graphene is a zero gap semiconductor. Wrap- 
ping the considered sheet of graphene, i.e., imposing pe- 
riodic boundary conditions (PBCs) around the circum- 
ference, yields a SWNT and leads to the formation of 
transverse subbands. For the low energy electronic struc- 
ture of metallic SWNTs, only the subbands touching at 
the Fermi-points are of relevance. In the following we 
consider armchair SWNTs of finite length and impose 
open boundary conditions (OBCs) at the two ends of the 
tube, i.e., that the wave function vanishes at the armchair 
edges. This condition mixes the two inequivalent Fermi 
points F = ±i^o from the underlying graphene first Bril- 
louin zone and yields the linear dispersion relation of the 
finite size SWNT shown in Fig. [2l It is characterized by 
two linear branches r = ± of slope ^hvp with the Fermi 
velocity vp ^ 8.1 • 10^ The allowed quasi-momentum 
values are given hy n = {n^ + A)7r/L, where e L 
is the tube length and A accounts for the fact that Kq 
may not be an integer multiple of tt/L. The kinetic part 
of the Hamiltonian, yielding the energy relative to the 



Ho = eo^r^ n^cl^^CraK + ^ rNra , (1) 

rcr ra 

where eo = Hvftt/L is the level spacing, and eA = eoA is 
the band offset energy. Finally c]^^^ creates an electron 
with momentum and spin a in branch r and the op- 
erator Nra counts the total electron number in branch r 
and of spin a. 

The interaction part of the Hamiltonian is given by 

aa' 

(2) 

where ^, ^"f" are fermion field operators and we use the 
Ohno-potential— , 

with Uo = 15 meV^^ and e lA — 2.4^ is the dielectric 
constant of graphene. In the next step we express the 3D 
electron operators in terms of the ID fermion-field^ 

V 2-L 

and obtain 

^^(r) = \^^sgn{F)ilJrFa{x)^ fpr(PpF{r) . (5) 

rF p 

Here F = ±i^o denotes the two independent Fermi- 
points, p = ± the two sublattices of graphene, and the 



3 



coefficients fpr of the sublattice wave function (PpF{r) 
are given by 1/ a/2 for p = + and — r/V^ for p = — . The 
sublattice wave function itself reads 



ReLc 



(6) 



where Nl is the number of graphene lattice sites iden- 
tified by the lattice vector and Lq denotes the 
graphene honeycomb lattice in real space. Furthermore, 
Xpz {r — R — Tp) is the Pz wavefunction of a carbon atom 
living on sublattice identified by the sublattice vector 
Tp. Upon integrating Eq. (|2]) over the coordinates radial 
to the tube axis, one eventually arrives at a ID interac- 
tion potential characterized by density-density and non 
density-density contribution^ so that the total Hamilto- 
nian reads 



Hq — Ho 



(7) 



With the help of bosonization^^ it is possible to diagonal- 
ize the density part Hq + Vpp. Eventually the bosonized 
and diagonalized Hamiltonian takes the forrn^: 



jSq>0 
J 



N-ra + (eo - u+) Nra + re^ 



• (8) 



Besides the ground state, it accounts for all the possi- 
ble fermionic and bosonic excitations of a SWNT. The 
bosonic excitations are described by the first term on 
the right hand side. The indices refer to total/relative 
{5 = +/— ) charge/spin {j = c/s) modes. The energies 
ejSq are given by 



jS c+ 



I, eo^g = c-,5+,s- 

with q = UgTr/L for Uq G Z and 



, (9) 



7o Jo 



X 4:C0s{qx) cos{qx') , (10) 



the contribution of the long-ranged density-density pro- 
cesses. Indeed C/'°"s(a;, x') = [C/'"*'^ + C/'°*<^^]/2 is the 
sum of the interaction potentials for electrons living in 
the same (intra) and different (sublattices): 



^intra/inter(^^ x') = J J d?r^d^r'^ 

X vlFW±pF'{r')v±pF'r)VpF{r)U{r- r') . (11) 



The second summand of (|8j) is the charging term with 
the charging energy Ec — VFg=o and also comes from the 
long range part of the Coulomb interaction. It counts 
the energy one has to spend to put = Xlrcr elec- 
trons on the dot, no matter what spin a G {t^i} oi" 
pseudospin r G {+, — } they have. The second fine of (j8j) 
starts with an exchange term favoring spin alignment. 
The exchange-splitting. 



J = 



Ed 



2N^ — 



-i2Ko{R^-K) 



X [U{R - R') - U{R -R' ^Tp- T-p)] , (12) 

being proportional to the difference of the Coulomb inter- 
action for electrons on the same and on different sublat- 
tices, accounts for the contribution of short range pro- 
cesses. The next term in ([8]) refiects the energy cost 
for adding electrons of the same spin band in the same 
branch, i.e., the Pauli-principle, where the correction u'^ 



-i2Ko{R^-K) 



R,R' 



X [U{R - R') + U{R -R' ^Tp- T_p)] . (13) 

Finally, the last term accounts for a possible band- 
mismatch, see Fig. [2j 

The eigenstates of Hq + Vpp are spanned by 



\N,m) 



n 

j5q 



rrijSq 



'^jSq- 



-\N,0) 



(14) 



Here N and rh denote the fermionic and the bosonic con- 
figuration, respectively, such that the state lA/", 0) has no 
bosonic excitation. The fermionic configuration is given 
by the number of electrons in each branch with a certain 
spin N — (N-^^ N-i^ N-^^^ N-^i). These eigenstates will 
be used to calculate the contribution of the non-density 
part of the interaction, i.e., {N^m\Vnpp\N' ^m'). Away 
from half-filling, they only couple states close in energy 
and one is allowed to work with a truncated eigenbasis 
(we check convergence of the results as the basis is en- 
larged). As shown by Yoshioka and Odintsov^^, for long 
SWNTs a Mott-insulating transition is expected to oc- 
cur at half-filling due to umklapp scattering. As found 
in Ref.^ umklapp processes acquire increasing weight as 
half-filling is approached also for finite size tubes, a pos- 
sible signature of the Mott instability, and the present 
theory breaks down. In recent experiments^^ the obser- 
vation of the Mott transition in SWNT quantum dots 
was claimed. 



B. Low energy spectrum away from half- filling 

The low energy regime is where the energies that can 
be transferred to the system by the bias voltage and the 



4 



.4t 



r=- r= + 



N^=4n+3 



N^=4n+1 



Nc=4n 



FIG. 3: Lowest lying states for fillings A^c = 4n + 1 and 
iVc = 4n + 3. For simplicity only the configuration of the last 
partially filled shell is shown. 



temperature stay below eo- This means no bosonic ex- 
citations are present, i.e., m = (0,0,0,0), and also no 
fermionic excitations are allowed, i.e., the four bands will 
be filled as equal as possible: \Nra — Nr'a' \ < 1 Vra, rV^ 
Our starting point are the eigenstates, Eq. ([II]), of the 
Hamiltonian in Eq. ([8]), which accounts for the kinetic 
and the density part of the full Hamiltonian. Now we 
have to split the examination into two cases. 
At first we consider states with total charge Nc equal to 
4n, 4n + 1 and 4n + 3. Those are unambiguously de- 
scribed by the fermionic configuration N because they 
are not mixed by the exchange effects. The only impact 
of the short-range interaction terms on these states is 
given by an energy penalty for double occupation of one 
branch r, a common shift for all eigenstates with fixed 
Nc G {4n, 4n + 1, 4n + 3}. Therefore we are left with^ 



Ef^ = \e,NI + min (iV^^, iV.J 



\^Nra -^iV-ra+ (eo -M+)iVra+reA 



(15) 



for the energy. If eA 7^ 0, states with the maxi- 
mum allowed number of electrons in the r = — branch 
will be the groundstates. For Nc = 4n the pseu- 
dospin branches r = ± are equally occupied, yield- 
ing an unique Nc — 4n groundstate. The corre- 
sponding configuration is taken as reference configura- 
tion for the Nc = 4n + 1, 4n + 2 and 4n + 3 cases. 
The lowest lying states for Nc G {4n + l, 4n + 3} 
are presented in Fig. [3l E.g., for the case Nc = 
4n + 1 we obtain four possible states corresponding 
to A/^ G {(n + 1, n, n, n), (n,n + l,n,n), (n, n,n + l,n), 
(n, n, n, n + 1)}. For simplicity we introduce for the 
states with an unpaired electron in the r = — branch 
the notation | tr)? I ir)- For electrons in the r = + 
branch we set I •, t), hi)- 

Analogously, neglecting exchange effects and setting 
eA = for the moment, the groundstates for the Nc = 
An -\- 2 filling are represented by the six states | t^t)? 
I ;,;), I t,;), | ;,t), | ti,-) and |-,t;), where, e.g., | t,t) 
means two electrons with spin ^ one on each branch 
— and +. Here the different fermionic configurations 
mix under the influence of the Vnpp processes and the 
groundstate structure will change dramatically due to 
off-diagonal contributions 
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TABLE I: The six lowest energy eigenstates for the filling 
Nc = 4n + 2 of an interacting SWNT. Due to short-ranged 
interactions there are three degenerate states of total spin 
S — h and three non- degenerate states of total spin S — 0. 



i,t) = 
■,n> = 



-J/2, 
J/2 . 



(16) 



Diagonalization of the interaction matrix yields the 
groundstate spectrum as it is shown in table HI The 
energies in the table are given relative to £^o,4n+2 = 



{2n^ + 2n + 1) (eo - ) - ^ {2n^ + 2n) + 2ix+n. It 



\EcNl 

is clear that the states 1^) and \h) will always be excited 
states, while the spin triplet, 5* = fi, is energy degener- 
ate. Now the question arises which states, the triplet or 
the I a) state, are the groundstate of the system. In accor- 
dance with table m the condition for a triplet groundstate 
is given by: 



(17) 



For a dielectric constant e = 1.4 it holds J = 0.72 A ^ 

a 

and = 0.22 A ^. Hence we find in terms of the level 
spacing eo and the tube diameter d\ 



< 0.4548 A ^. 

a 



(18) 



Obviously this makes the triplet groundstate more un- 
likely compared to the S = ^ groundstate as it can be 
seen in Fig. |4j For a (6,6) nanotube of 300nm length, 
the band-mismatch must be eA < 0.3 meV = 0.06eo to 
be in a triplet groundstate. In the experiments^^^^ band- 
mismatches are of the order of 0.3eo and, as expected 
from our theory, \a) - groundstates are observed. 



III. SPIN-DEPENDENT TRANSPORT 

In this section we discuss the set-up to evaluate spin- 
dependent transport across a SWNT weakly coupled to 
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FIG. 4: Phase diagram to determine the groundstate of dif- 
ferent tubes of length 300nm. The chance to find a triplet 
groundstate increases with increasing exchange parameter J, 
i.e., with decreasing tube diameters. 



leads, see Fig. [TJ and the main calculation tools. The 
Hamiltonian of the full system reads 



E 



Hi+Ht 



Hp. 



(19) 



where I = s^d denotes the Hamiltonian in the source and 
the drain contact, respectively. The leads magnetization 
is accounted for in terms of a Stoner Hamiltonian where 
the density of states, Vi^{e)^ for the majority (a =t) 
and the minority (a =\) carriers are different. We treat 
the leads within the wide-band approximation, i.e., we 
regard the density of states as constant quantities to be 
evaluated at the leads chemical potentials (is and fid- We 
can thus define the polarization by (/ = s^d): 



(20) 



Moreover, we will consider a symmetric set up Vsa = 
Vda = ^cr and Ps = Pd = P' The total density of states 
is given by Vtot = + account for the bias volt- 

age Vb in terms of the difference eVb = jis — IJ^d between 
the electrochemical potentials in the source and drain 
leads. Further, Ht in Eq. (p!9|) is the tunneling Hamil- 
tonian which we will treat as a perturbation since weak 
coupling to the leads is assumed. Finally, He^t describes 
the influence of the externally applied gate voltage Vg. 
The gate is capacitively coupled to the SWNT and hence 
contributes via a term eaVgNc with a a proportionality 
factor. 

In order to evaluate the current-voltage characteristics 
we use the method developed in Ref.^^ where, starting 
from the Liouville equation for the density matrix of the 
full system, a generalized master equation (GME) for the 
reduced density matrix p (RDM) of the SWNT is ob- 
tained to second order in Ht- Once the stationary RDM 



parameters 


label 


value 


length 


L 


300.06 nm 


diameter 


d 


0.81 nm 


dielectric constant 


e 


1.4 








charging energy 


Ec 


6.7 meV 


level spacing 




5.6 meV 


Coulomb excess energy 




0.15 meV 


exchange energy 


J 


0.49 meV 


orbital mismatch 




OmeV or 1.68 meV 


thermal energy 




4.0 X 10"^ meV 


transmission coefficient 




1 X 10"^ meV 



TABLE II: Parameter set of the 300nm (6,6) nanotube inves- 
tigated in this work. 



is known, the stationary current through e.g. the source 
lead is evaluated from the relation Is = eTr{pNs}^ where 
Ns is the number operator for electrons in the left lead. 
As this procedure with the relevant equations is thor- 
oughly explained in Ref.^^, we refrain from repeating it 
here. The GME can be solved in analytic form in the 
linear regime, being the focus of the following Sec. IV. 
In the nonlinear regime, discussed in Sec. V, the differ- 
ential conductance is evaluated numerically. Moreover, 
from here on we will focus on the transition between 
charge states 4n + 1 < — > 4n + 2, mirror symmetric to 
4n + 2 i — > 4n + 3, as these two transitions are the ones 
that reveal exchange effects. The remaining transitions 
4n i — > 4n + 1 and 4n + 3 i — > 4(n + 1) will not qualita- 
tively change due to the presence of short range processes 
and we hence refer to the discussion in^^. 
If not otherwise specified, we choose nanotubes described 
by the parameters in table [III In order to obtain an \a) 
groundstate we assume a band- mismatch of eA = 0.3eo = 
1.68 me V, whereas for a triplet groundstate we choose 
eA = 0. 



IV. THE LINEAR REGIME 
A. Conductance at zero magnetic field 

We focus on the conductance formulas for the two cases 
of tunneling from the 4n + 1 groundstates into the S = 
groundstate \a) or into the triplet groundstates. 

For the transition |cr, •) < — > \a) the conductances in the 
case of parallel, G^, and antiparallel, G^^, magnetized 
leads are found to be 



2 2 



7 



G: 



AP 



n 1 + 7 

(P2_ 1)^(1 + ^) 



2 - /(/ia) 



2 -> 



(21a) 
(21b) 



a) 



b) 



a) 



b) 



P=0 
P=0.4 
P=0.8 
P=0.99 




331.94 331.96 331.98 332.00 332.02 

Vg [meV] 



Source 



Drain 



c) 



d) 




331.94 331.96 331.98 332.00 332.02 

Vg [meV] 



Source 



Drain 



FIG. 5: Panels a) and c). Conductance vs. gate voltage for 
the |cr, •) i — > I a) resonance for parallel, , and antiparal- 
lel, Ga^^ lead magnetization. In both cases the analytical 
predictions Eqs. (|21a|) . (|21b|) (continuous curves) perfectly 
match with the results from a numerical evaluation of the 
GME (squares). Strikingly Ga is independent of the leads 
polarization P, while Ga^ is maximal at P = 0. Panels 
b) and d). Schematic explanation of the different polariza- 
tion dependence. The red spin specifies the spin of the state 
|cr, •). The dashed/continuous arrows indicate rare/favorable 
tunneling processes. For parallel magnetization, panel b), the 
fast tunneling channel is the one with an excess spin ^ and the 
electron transferred from source to drain is always a majority 
electron t- If the initial dot spin is this is likely to tunnel 
to the drain, such that at the end of the tunneling process a 
spin- flip has occurred, leaving the dot in the favorable con- 
figuration with a spin For antiparallel lead magnetization, 
panel d), the fast channel corresponds to one electron in the 
dot with spin t- To this channel, however, is associated a spin 
flip. Because the situation with initial spin \. involves a rare 
tunneling process from the source lead, the conductance gets 
diminished by increasing polarization. 



with c = ci/y^Ci + 1, the Fermi function /(/i) eval- 
uated at the gate voltage dependent energy difference 
jiia = £^|a) ~ ^|cr, ) and P the inverse temperature. The 
parameters <l> = <l>5 and j = ^d/^s describe the possible 
asymmetric lead transparencies^^ (hereby, $ is in second 
order of the tunneling coupling contained in Ht)- The 
conductances are shown in Fig. [5^) and[5fc) for the sym- 
metric transparencies case 7 = 1 and Vtot^ = 10~^meV. 
Strikingly, in the parallel magnetized case there is no de- 
pendence on the polarization since there is never a block- 
ing state involved in transport, see Fig. [SJd). For the 
antiparallel case, in contrast, transport is limited by the 
weakest channel (when there is a | - electron on the dot) 
and one can drive the conductance to zero by tuning the 
polarization to P ^ 1. This feature is explained in Fig. 
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FIG. 6: Panels a) and c). Conductance vs. gate voltage 
at zero band-mismatch (triplet groundstate) for parallel, , 
and antiparallel, Gf^, lead magnetization. G^ is indepen- 
dent of the leads polarization P, while Gf^ is maximal at 
P = 0. The absolute value of the conductance is slightly 
larger than for the |cr, •) i — > \a) case since more channels are 
involved. Panels b) and d). Schematic explanation of the dif- 
ferent polarization dependence. For simplicity we only drew 
the case in which the initial excess spin (red spin) is in the 
r — -\- branch. For parallel magnetization, panel b), the fast 
channel corresponds to the It?*) < — > transition which 

conserves the spin of the excess dot electron. For antiparallel 
magnetization, panel d), the fast channel corresponds to an 
initial excess spin t electron likely to tunnel to the drain and 
being replaced by a spin ^ from the source. The situation with 
an initial spin ^, however, corresponds to a weak channel. In- 
creasing the polarization highly populates the |t-i) state and 
transport decreases. 



For the case of the triplet groundstate we face a com- 
pletely new situation. First, we have for A^c = 4n + 1 
filling four degenerate states available because the band- 
mismatch has been chosen to be zero. Secondly, we cou- 
ple to three different states in the case of Nc = 4n + 2 
rather than to just one. However, the conductance plots 
do not qualitatively change as it may be seen in Fig. [6^1) 
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andEt). The conductance formulas read: 



a) 



r^AP 



n 1 

(p2_ 



P Vtot ^ 



-7 
1)7(1 



f{/^t)f{-fJ^t) 



-7) 



4 -/(/it) 



(22a) 
(22b) 



Compared to Eqs. (j2Taj), (|21b[) the prefactor changed 
from to 3 due to the three involved triplet states. 
The quantity /it = E\t^ — Ei is the difference between 
the triplet and the A/'c = 4n + 1 - groundstate energies. 
In addition, the denominator in the term containing the 
Fermi-functions has also changed to account for the de- 
generacy of the 4n + 1 - filling states. The qualitative be- 
havior, however, does not change compared to the case of 
an I a) groundstate, such that one cannot determine the 
spin nature of the groundstate from these plots alone. 



B. Conductance in the presence of an external 
magnetic field 

In this section we consider the infiuence of an exter- 
nally applied magnetic field (Zeeman-field) which clearly 
reveals the character of the groundstate for 4n + 2 and, 
moreover, may even change the groundstate depending 
on the field strength. The field causes an additional Zee- 
man energy to states with a spin-component Sz ^ 0. 
The sign is negative if the concerned state in the tube 
is parallel to the external field and positive if antiparal- 
lel. Thus, the chemical potential differences appearing 
in Eqs. (|2Ta|) . (|2Tbl) . (|22a|) and (j22bl) wih be shifted by 
±Ez = ±/iB^- We use the convention ji^ = — Ez and 
/i| = ji^Ez. Furthermore, in order to improve the read- 
ability, we introduce the abbreviation f±^/^ = /(±/i>^/|). 
The conductances for the antiparallel set-up are 
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We do not find qualitative differences with respect to the 
zero magnetic field case: the conductances decrease in 
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FIG. 7: a) Conductance near the |cr, •) < — > |a) transition for 
parallel magnetized leads and applied magnetic field. The 
peaks corresponding to higher polarizations are shifted to 
lower gate voltages, b) Schematic explanation of the polariza- 
tion and gate-voltage dependence for small (left sketch) and 
large (right sketch) polarization. The red spin indicates the 
spin of the excess electron initially present on the dot. The 
thick and thin lines are frequent and less frequent transitions, 
while dashed lines indicate rare transitions. Large polariza- 
tions favor processes involving majority spins while, due to 
the extra required Zeeman energy, the Fermi function sup- 
presses processes where a spin ^ is transferred. Thus at small 
polarizations the transport is mostly mediated by spin \, - elec- 
trons while at large polarizations t - electrons are preferred. 
Correspondingly the peak position is shifted to smaller gate 
voltages as the polarization is increased. 



both cases with increasing polarization. In the following, 
we will therefore only focus on the parallel case, where 
we find interesting behavior for small Zemann splittings. 
The conductance formulas for parallel lead magnetization 
take the form 
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The corresponding plots can be seen in Figs. and[8K). 
In these calculations we considered a small magnetic field 
of 0.07 T which equals in magnitude the thermal energy of 
ksT = 0.004 meV. This provides a situation with a finite 
occupation probability for all included states. Specifi- 
cally, this means that also states containing ^ - electrons 
will be populated, but the population of states contain- 
ing t - electrons will be preferred. The first thing we 
observe in both Fig. [7^) and [8^) is that the once degen- 
erate curves in Figs. [5^) and [6^) now split into distinct 
curves for the four different polarizations. Moreover, the 
peaks of the curves corresponding to less polarized leads 
continuously move to higher gate voltages. Finally the 
conductance decreases /increases with increasing polar- 
ization for the a/t cases, respectively. Let us examine 
the results starting with the \a) - groundstate. We will 
divide the analysis in two cases, slightly polarized leads 
and strongly polarized leads. 

For only slightly polarized or non-polarized leads the sit- 
uation is intricate as we have to deal with competing pro- 
cesses. On the one hand there is a highly populated 
I ^7 •) state and a slightly populated | i, •) state in the 
tube. From this point of view, the system prefers | - 
electrons to tunnel into the \a) state and to leave the 
dot subsequently such that the tube always remains in 
the preferred | ') state (Fig. [71 sketch b), upper left 
panel). Only rarely, the t - electron tunnels out, as this 
would result in a spin- flip to the disfavored | i, •) state 
(Fig. [71 sketch b), lower left panel). On the other hand, 
entering of | - electrons is suppressed compared to trans- 
port of t - electrons, not so much by the small polar- 
ization, but mainly due to the Zeeman splitting in the 
involved Fermi- functions: The chemical potential for ^ - 
electrons exceeds the one for t - electrons by 2Ez such 
that > at any gate voltage. However, in the end 
it will be a mixture of mainly ^ - electrons and some t 
- electrons responsible for transport. This can also be 
seen by the fact that the curves for small polarizations 
are shifted to higher gate voltages which accounts for the 
higher chemical potential of the ^ - electrons. In addition, 
the total amplitude of the conductance is decreased com- 
pared to the case without the magnetic field. Fig. [S^i), 
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FIG. 8: a) Conductance near the triplet resonance for parallel 
magnetized leads and applied magnetic field. In contrast to 
the case of a singlet resonance, Fig. [71 transport increases as 
the polarization is enhanced, b) Schematic explanation. At 
small leads polarization the distribution of t - electrons and 
^ - electrons is almost equal. However, the |ti) - channel is 
preferred to the others. Increasing the polarization enhances 
the dominance of this channel and correspondingly the con- 
ductance. Simultaneously the conductance peak is shifted to 
lower gate voltage indicating the dominance of t - electrons. 



as there is always a limiting element - either the small 
Fermi-function or the small population - involved. 
In the case of highly polarized leads we face the situation 
where there are very few ^ - electrons in the leads. As 
temperature provides a small, but nonzero population of 
the slightly excited state | |, •), current mainly flows via 
the polarization-favored t - electron channel. Since the 
chemical potential, the increment of the Fermi- functions, 
is smaller than in the former case the transition takes 
place at slightly lower gate voltages. The situation again 
is visualized in the sketch b) of Fig. [71 in the upper and 
lower right panel. 

At the triplet resonance we observe not only quantita- 
tive, but also qualitative changes. The plot can be seen in 
Fig. [H^i) and all relevant tunneling processes are sketched 
in Fig. [Hb). Let us again start with unpolarized or just 
slightly polarized leads. Due to a large population of the 
spin t states in the Nc = 4n + 1 case and of the \ti) state 
in the A^c = 4n + 2 case transport is mainly mediated 
via the majority charge carriers, i.e. t - electrons (Fig. 
[8)3), upper right panel). However, the resulting current 
is smaller than in the case without magnetic field since 
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it is harder to make use of the I - electrons that are stih 
largely at disposal in the leads. 

A high polarization decreases the number of ^ - elec- 
trons in the leads in favor of the t - electron number, 
and such transport via the already preferred \ti) chan- 
nel is strongly enhanced. As a consequence, the conduc- 
tance by far exceeds the conductance without magnetic 
field and polarization. This effect should be detectable 
in an experimental setup and would give a possibility to 
distinguish between a triplet groundstate and a 6* = 
groundstate. 



V. THE NONLINEAR REGIME 

In the finite bias regime also excited states become 
available and, due to the resulting high number of in- 
volved states, it is necessary to calculate the current 
numerically. We show the current and the stability di- 
agrams - the differential conductance '^^(Vb^Vg) as a 

function of the gate and the bias voltage. The stabil- 
ity diagrams give a clear indication whether the involved 
groundstate in the transition 4n + 1 < — > 4n + 2 is the 
I a) state or the triplet. In the case of antiparallel lead 
magnetization we find negative differential conductance 
(NDC) for transitions involving the \a) state. We also 
observe NDC for transitions involving the \a) state or 
the triplet if an external magnetic field is applied. 
The current as a function of the gate and the bias voltage 
is shown in Fig. [9ti) for the \a) groundstate and in Fig. 
[9)3) for the triplet groundstate. All states with up to 
one bosonic excitation have been included in the calcula- 
tion. A 4-electron periodicity of the Coulomb diamonds 
is clearly seen. The change in color indicates a change in 
current and therefore the opening of a new channel. At 
high bias a smearing of the transitions due to the multi- 
tude of bosonic excitations is observed. In the remaining 
of this section we focus on the gate voltage region rele- 
vant for the 4n + 1 < — > 4n + 2 transitions. In the plots of 
the differential conductance reported in the following we 
did not include the bosonic excitations to avoid a multi- 
tude of transition not relevant for the coming discussion. 
A polarization P = 0.9 is chosen. 



A. Differential conductance at zero magnetic field 

Figs. [TOk ) and Fig. [TOb) show the stability diagrams 
for parallel and antiparallel lead magnetization, respec- 
tively, for the case of the \a) groundstate. The two tran- 
sition lines h and e were emphasized by a dashed line 
because these lines are so weak that it was not possible 
to resolve them together with the other stronger lines. 
The most obvious difference between the parallel and 
the antiparallel setup is the weakness of all transition 
lines beyond the triplet occupation (line b) for antipar- 
allel lead magnetization. Moreover an NDC line, (line 
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FIG. 9: Current versus gate and bias voltages for unpolarized 
leads. In total 176 states have been included, which corre- 
sponds to all states with at most one bosonic excitation. For 
4n + 2-filling this amounts to 32 different states, a) Band- 
mismatch eA = 0.3 eo corresponding to an 5* = groundstate 
for the 4n + 2 filling, b) Band-mismatch ca — correspond- 
ing to an S = h groundstate at filling 4n + 2. In both cases a 
4-electron periodicity of the Coulomb diamonds is observed. 



6), not present in the parallel magnetization case, is ob- 
served. 

In order to explain the line positions in Fig. [TQk).b) 
we provide a schematic drawing in Fig. [11] which is 
based on a bias trace at the particular gate voltage which 
aligns the groundstates (white vertical lines in Fig. [TQ|) . 
The differently colored arrows stand for new transport 
channels that open at certain bias voltages. The chan- 
nels open in the order of a to e for transitions from 
4n + 1 — > 4n + 2 (dashed arrows) and f to h for transi- 
tions from 4n + 2 — > 4n + 1 (solid arrows). Sometimes 
opening of a new channel also opens other channels that 
have been blocked before and one does not see distinct 
lines for these transitions. Fig. [11] relates the concerned 
transitions to the required bias voltages. Moreover, the 
line g stands for transitions between the triplet and the 
I -,(7) states, i.e., it is a transition between excited states. 

To explain the NDC in Fig. [TOb) which follows upon 
line b in the range between lines / and line ^, we ob- 
serve that -in correspondence of the b line- below the 
resonance only the transitions from •) to the \a) state 
is possible. Above resonance also the triplet \t) is ac- 
cessible. For the case of antiparallel polarization, both 
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FIG. 10: Differential conductance for transitions between 
4n + 1 i — > 4n + 2 filling in the \a) - groundstate. The po- 
larization has been chosen to be P — 0.9. The four lowest 
lying states for 4n + 1 and the six ones for 4n + 2 filling were 
included. The vertical white line is the bias trace we follow 
to explain the distinct transition lines in Fig. 1111 a) The 
leads are magnetized in parallel, b) Antiparallel magnetized 
leads. We observe a different intensity of the excitation lines 
between parallel and antiparallel magnetization. In partic- 
ular a pronounced negative differential conductance (NDC) 
occurs in correspondence of the transition between |cr, •) and 
the triplet (line h). 
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FIG. 11: Schematic drawing for the possible transitions oc- 
curring by sweeping the bias voltage at the gate voltage that 
aligns the |cr, •) and the | a) -states (white dashed line in Fig. 



filling groundstates. Line a is the groundstate to ground- 
state transition. Lines b to d indicate transitions from 
the 4n + 1 groundstates to |a), \s) and respectively. 
They come in the expected order, at an applied voltage 
Vb/2 equal to J and J + as it is shown in table 
[H Line e stands for the transition from the triplet to one 
of the 4n + 1 groundstates. 

For the antiparallel setup, Fig. [T2b), we may see the 
same effect as we have observed in Fig. [TOb). i.e., all 
lines beyond the transition to the triplet decrease in in- 
tensity. Since the triplet is the groundstate, this means 
all excitation lines are weak and may not be resolved in 
the figure. 



B. Differential conductance in parallel magnetic 
field 



provide only weak transport channels: below the reso- 
nance transport is mostly mediated by t - electrons (see 
also sketch of Fig. [5j) which are minority electrons for 
the source contact; above resonance, after some tunnel- 
ing processes the system will always end up in the \t-i) 
state which is a trapping state. Just at the exact reso- 
nance, the thermal energy allows electrons to tunnel forth 
and back, i.e., a | - electron has the possibility to tun- 
nel back into the source contact and transport is slightly 
enhanced. Once the bias voltage exceeds the exact res- 
onance the trapping state gets occupied for long 
times and the current diminishes again. 

The fact that the | i, •) < — > transition serves as 

the major transport channel once it has been opened is 
also the reason why all transition lines above line b are 
so weak. 

In Figs. [T2b ) and [T2b) the stability diagrams for the 
S = H triplet groundstate are shown. They look a lot 
simpler than the ones in Fig. [TOldue to the absence of a 
band-mismatch, causing a degeneracy of all four 4n + 1 



Here we present results for an applied magnetic field 
of strength =0.1 meV, Fig. [131 The leads are par- 
allel magnetized and a polarization of P = 0.6 has been 
applied. The magnetic field removes the spin degeneracy 
of the triplet as well as of the 4n + 1 filled states; the 
resulting Zeeman split transitions are clearly seen in Fig. 
[13] a) and are less well resolved in Fig. [13] b). 

Explicitly, for the \a) groundstate, line b from Fig. [10] 
splits into lines b and c in Fig. [13] We notice that line c 
shows an NDC effect due to the opening of the channel 
II,-) — > though this transition, as mediated by 

minoriy | - electrons, is rare, once it happens the system 
is trapped in the \t-i) state for a long time due to the 
parallel polarization of the leads. For transitions from 
4n + 2 to 4n + 1 line A: is a new line that was Coulomb 
blocked in Fig. [10] It denotes the transition 1^) — > |-,|) 
and ends in line e since the \s) state must be populated. 
Also, we notice the absence of the \s) — > liiie since 
it is Coulomb blocked by the groundstate to groundstate 
transition (line j). 

For the S = ft triplet groundstate. Fig. [T3b). we ob- 
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FIG. 12: Differential conductance for transitions between 
4n + 1 i — > 4n + 2 filling in the triplet groundstate. The 
polarization has been chosen to be P — 0.9. The four lowest 
lying states were included for 4n + 1 and the six lowest ones 
for 4n + 2. a) Leads parallel magnetized, b) Leads polarized 
antlparallel. From the stability diagrams it is possible to di- 
rectly extract the exchange parameters u'^ and J since the 
bias voltage 14/2 = u'^ is needed to open transition line b 
and V6/2 = J to open line c. 



serve that line b and line c show NDC effects. Line b rep- 
resents transitions from — > \to) or | •) — > |to), 
which is not a trapping state. However, the applied bias 
voltage is sufficient to also populate the |-,|) and | i, •) 
states from \to) and subsequently from |-,^) and | |, •) 
the trapping state This process is also visualized 

in Fig. [lU In the very same way it is possible to get 
trapped in the state via the \a) state indicated by 

line c. 



C. The magnetic field sweep 

In a seminal experiment Moriyama et al.^ demon- 
strated a transition from a 6* = groundstate to di Sz = h 
groundstate upon magnetic field sweep in a SWNT quan- 
tum dot. In this section we have computed the differen- 
tial conductance in a gate- voltage and magnetic field plot 
both for unpolarized, as in^, and parallel polarized leads 
with P 0.9. 

We start from the \a) groundstate at 5 = with a 
band- mismatch of 0.24 eo (smaller than we previously 
used). This choice yields a change of groundstate from 
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FIG. 13: Differential conductance for transitions between 4n+ 
1 < — > 4n + 2 filling with an applied magnetic field of = 0.1 
meV. A parallel lead magnetization was assumed with the 
polarization P — 0.6. a) |a) - groundstate. Soon after line c 
an NDC effect is observed due to the occupation of the \t-\) 
trapping state, b) Triplet groundstate. After lines b and c 
NDC occurs due to an increased population of the \t-\) state. 



I a) to the triplet at a magnetic field 6 T as measured 
experimentally^. To observe well visible patterns, we in- 
creased the temperature by a factor of ten compared to 
Tab. [m 

The result of our calculation is presented in Fig. [T5k). 
At a gate voltage of approximately 0.322 meV and 
0.323 meV we have two V-shaped transition patterns (a 
and 6) each of width 2Ez = 2\ibB. The separation be- 
tween a and h at zero field is the band- mismatch eA- In- 
terestingly, for polarized leads, the branches belonging to 
transitions involving (| i, hi))^ corresponding to the 
positive slope of the "V", are NDC lines. Fig. [T5b). The 
reason is the same as addressed already in section IV Bl 
once the ^ - channel becomes available, there is some 
chance that from time to time a minority charge car- 
rier {\ - electron) enters from the source. As the drain 
is polarized in parallel to the source, it will take quite 
a while until this electron can leave the SWNT again, 
such that transport gets hindered. At the gate voltage 
of approximately 0.328 meV, one enters the A^c = 4n + 1 
Coulomb diamond (line c) and transport gets completely 
suppressed. The dot is in the groundstate | t, •) at 5 ^ 0. 
At Vg 0.329 meV transport from A^c = 4n + 1 to the 
I a) state is enabled (line d). 

The next transitions (patterns g) we observe are 
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FIG. 14: Schematic drawing of the possible transitions if the 
|-,t) and I t, •) states are ahgned to the |ti) state by the gate 
voltage at finite magnetic field and in the triplet groundstate. 
It provides the explanation for the transitions lines observed 
in the inset of Fig. [T3bV 



again split by 2Ez and therefore shaped like a "F" . In 
all cases, the positively sloped branches are now again 
of NDC nature for a parallel lead polarization. The first 
"V" belongs to the triplet (pattern e) and is of stronger 
intensity than the following two patterns. The transitions 
I •) ^ — ^ 1^1) cind I ^, •) < — > |to) contribute to the 
negative sloped part, while | ') ^ — ^ l^o) and | i, •) < — > 
are responsible for the positive shaped line. The 
crossing of the e and d lines occurring at 5 = 6 T, point 
P, indicates the change in the groundstate from \a) to the 
state \ti). 

From the triplet pattern e the additional gate voltage 
equal to the exchange energy J is needed to arrive at the 
last two "F" - shaped patterns / and g. Compared to 
the lines for the triplet transition they are quite close to 
each other and of less intensity. These lines belong to a 
transition from both the | i, •) and the | ') states to 
the \s) - singlet (pattern /) and the \b) state (pattern 
g). Finally, the lines on the right edges of the plots are 
mirror images and belong to backward transitions from 
A/'c = 4n + 2 to A/'c = 4n + 1; for this reason they mark 
a decrease of current for both polarized and unpolarized 
leads. 
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FIG. 15: a) Differential conductance dl/dVg for a 5-field 
sweep in the \a) - groundstate case. The applied bias voltage 
was fixed at 5.8 meV. Red lines indicate transitions that be- 
come possible at a certain gate voltage and blue lines show a 
transition that drops out of the transport window. The "V"- 
shaped patterns a and b represent transitions from Nc = 4n 
to |cr, •) and |-,cr), respectively. Each of the patterns is split 
by 2Ez denoting t - electrons and ^ - electrons tunneling in. 
At line c we enter the Nc = 4n + 1 Coulomb diamond and 
transport gets suppressed. Line d stands for the groundstate 
to groundstate transition from | •) to the \a) - state. The 
"V" -shaped pattern e is due to the transition Nc = 4n + 1 
to the triplet whereas / and g denote transitions to the \s) 
- singlet and the \b) state, respectively. At the point P the 
groundstate changes from the \a) state to the |ti) - triplet, b) 
Ferromagnetic leads, polarized in parallel with P = 0.9, are 
assumed. This changes the intensity of the transitions, while 
their positions are preserved. Moreover, transitions to excited 
states involving spin-down electrons are disfavored channels 
and hence converted from positive to negative differential con- 
ductance lines. 



VI. CONCLUSIONS 

In summary, we have calculated spin dependent trans- 
port through fully interacting SWNTs in both the linear 
and the nonlinear regime, with and without an applied 
magnetic field. 

Peculiar of metallic SWNTs of small diameter is the 
possibility, due to exchange interactions, to find the sys- 
tem at 4n + 2 filling either in a groundstate of total spin 
S = or S = h. Which of the two groundstates occurs 
in a real nanotube depends on the relation between the 
exchange energy and the orbital band mismatch. Thus, 



with focus on transitions involving 4n + 1 i — > 4n + 2 fill- 
ing, we investigated both situtations and demonstrated 
pronounced differences in the current-voltage character- 
istics depending on the considered groundstate. 
For example in the linear regime the conductance for 
parallel lead magnetization and finite magnetic field in- 
creases by raising the polarization for the case of a triplet 
groundstate but it decreases for the S = groundstate. 
This is due to the fact that for the triplet groundstate 
transport is dominated by a channel involving the triplet 
state \ti) (with both spins t); for the 5* = case trans- 
port to be mediated by the majority electrons requires 
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to make use of the 4n + 1 lowest excited state | •) (and 
hence less favorable), Zeeman split from the ground state. 

In the nonlinear regime we presented stability dia- 
grams with parallel and antiparallel lead magnetization 
for both ground sates. In the antiparallel case it was 
possible to observe a negative differential conductance 
(NDC) effect for the S = groundstate, following imme- 
diately upon a conductance enhancement at the opening 
of a trapping channel to the excited triplet state 
Directly at that resonance, electrons can, just by thermal 
activation, tunnel back and fourth, such that trapping in 
the \t-i) state can not yet act, leading to an intermediate 
conductance increase. Away from resonance, the block- 
ing effect fully occurs, resulting in the NDC. By adding 
an external magnetic field in the parallel setup we found 
NDC effects for both groundstates caused by spin block- 
ing mediated by | - channels, involving in particular the 
triplet state 

Finally, we also presented results for the differential 
conductance in a gate-voltage and magnetic field map 
at finite bias. These magnetic field sweeps immediately 
allow to recognize the nature of the 4n + 2-filling ground- 
state at zero field, as well as to tune the nature of the 
groundstate from S = to Sz = h upon variation in the 



field amplitude. Our results for unpolarized leads are 
in quantitative agreement with experiments on a small- 
diameter SWNT by Moryama et al.^. Importantly the 
sweep at zero field also allows to immediately read off the 
values of the short range interactions J and . Specif- 
ically, J is the singlet-triplet exchange splitting and 
characterizes at zero orbital mismatch the energy differ- 
ence between two of the low energy states of total spin 

= 0. In the presence of polarized leads the magnetic 
field sweep also reveals lines of NDC due to the trapping 
nature of all ^ - channels. 

The predictions of our theory are in quantitative agree- 
ment with experimental results obtained so far for unpo- 
larized leads^— . Due to recent achievements on spin- 
polarized transport in SWNTs^^^^, our predictions on 
spin-dependent transport are within the reach of present 
experiments. 
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